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Abstract 



In this article, we show a result of approximation in law to subfractional Brownian 
motion, with H > h, in the Skorohod topology. The construction of these approxi- 
mations is based on a sequence of I.I.D random variables. 
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1. Introduction 

Recently, the long-range dependence property has become an important aspect of stochas- 
tic models in various scientific areas, such as hydrology, telecommunication, finance and 
so on. The best known and widely used process that exhibits the long-range dependence 
property is the fractional Brownian motion (FBM) introduced by Mandelbrot and Van 
Ness [9]. FBM is a suitable generalization of the standard Brownian motion. It exhibits 
long-range dependence and self-similarity. Moreover, it has stationary increments. Refer 
to Samorodnitsky and Taqqu [12] for more information on FBM. 

On the other hand, many authors have proposed to use more general self-similar Gaus- 
sian processes and random fields as stochastic models. Such models have raised many 
interesting theoretical questions. However, contrast to the extensive studies on FBM, 
there has been little systematic investigation on other self-similar processes. It seems 
to me that the main reason for this is the complexity of dependence structures for self- 
similar processes, which do not have stationary increments. Therefore, it seems interesting 
to study these processes. 

Recently, Bojdecki et al. [3] introduced and studied a rather special class of self-similar 
processes which preserves many properties of FBM. This process arises from occupation 
time fluctuations of branching particle systems with Poisson initial condition. This process 
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is called the sub- fractional Brownian motion (Sub-FBM). It is also independently appeared 
in Dzhaparidze and Van Zanten [5]. 

Weak convergence to the FBM processes has been studied extensively since the works of 
Davydov [4] and Taqqu [TJ]. In recent years, many new results on approximations of FBM 
have been established. See, for example, Meyer, Sellan and Taqqu [10], Li and Dai [8] and 
the references therein. When we consider the Sub-FBM process, a natural question is to 
ask what processes can approximate to it. In recent years, some results on approximations 
of Sub-FBMs have been obtained. For example, Bardina and Bascompte pQ presented a 
weak theorem for Sub-FBM based on a Poisson process. Harnett and Nualart [7] proved 
weak convergence to Sub-FBM with H = ^. Garzon et al. [6] proved a strong uniform 
approximation with a rate of convergence for Sub-FBM by means of transport processes. 

Inspired by these results, we will present a approximation to Sub-FBM with H > ^. 
The paper is organized as follows. In Section 2, we recall some preliminaries for Sub-FBM 
and present the main result of this paper. Section 3 provides the proof of the main result. 

2. Preliminaries and Main result 

In this section, we first briefly recall some basic preliminaries for Sub-FBM. The sub- 
fractional Brownian motion X H = {X H (t),t > 0} with index H G (0, 1) is a centered 
Gaussian process with covariance function 

E[X H (t)X H (s)] =s 2H + t 2H - l -[{s + t) 2H + \t-s\ 2H ]. (2.1) 
Furthermore, it satisfies the following estimates for s < t 

[(2 - 2 2H ' X ) A l](t- s) 2H < E[(X H (t) - X H (s)) 2 ] < [(2 - 2 2 "" 1 ) V l] 2H (t - s) 2H (2.2) 

Sub-FBM X H is neither a semimartingale nor a Markov process unless H = ^. When 
H = i, we have the standard Brownian motion W. Sub-FBM has properties analogous 
to those of FBM (self-similarity, long-range dependence and Holder paths). But its incre- 
ments are not stationary. More works on Sub-FBM can be found in Tudor [15], Yan and 
Shen \17\ [T8] and the references therein. 

In this paper, we assume the index H G (i, 1). Using the Hankel transform in Dzha- 
paridze and Van Zanten [5], we get from Tudor [16] that the process 

W(t)= f ^ H (t,s)dX H (s) (2.3) 
Jo 

is the unique Brownian motion such that 

X H (t) = I K H {t,s)W(ds), t e [0, 1], (2.4) 
J o 

where 



<Mi, s) = -^4- \t H -l (t 2 - s 2 ) ^ H ~(H-h f (x 2 - s 2 ) *- H x H -idx] 1 (0 , t) (s), 
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and 

with Ch being the normalizing constant. Moreover, X H and W generate the same filtra- 
tion. 

Let us recall some known facts. Consider a sequence of I.I.D random variables 
with E[&] = and E[^ 2 ] = 1. The Donkser invariance principle states that the sequence 
of the processes 

[nt\ 

W n (t)=J2tk (2.5) 
k=i 

converges weakly to a standard Brownian motion in the Skorohod topology. Here [^J 
stands for the greatest integer not exceeding x. 
This result has been extended by Sottinen [13] to the fractional Brownian motion. Define 



FJt,s)=n( F(^-,u)du, n>l (2.6) 
J,-± n 

n 

where F is the kernel that transforms the standard Brownian motion into a fractional 
Brownian one, i.e., 

1 1 /"* 1 3 
F(t,s) = c H (H- / u H ~^(u- s) H -^du, 

2 J s 

where ch is still the normalizing constant. 
Set 

Z n {t)= / F n (t,s)W n (ds) = Vn /" F(^,u)du-^=. (2.7) 
Jo Ji=± n V n 

Then {Z n (t)} converges weakly to FBM (see Sottinen [T3] ). 
Inspired by Sottinen [13], we define: 

K n (t,s)=nf K H (^-,u)du, n > 1. (2.8) 

n 

Then K n (t, •) is an approximation of Kn(t, •) for every t G [0, 1] 
Set 

X n (t) = n K n (t,u)W n (du) = Y y n " K^(i^, M )dn^. (2.9) 

In this paper, we will prove 

Theorem 2.1 T/ie /aws of processes {X n (t), t £ [0, 1]} converge weakly, as n tends to 
infinity, to the law of the sub -fractional Brownian motion X H given by {2.1$. 



In the rest of this paper, most of estimates contain unspecified constants. An unspec- 
ified positive and finite constant will be denoted by C, which may not be the same in 
each occurrence. Sometimes we shall emphasize the dependence of these constants upon 
parameters. 
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3. Proof of Theorem [2J] 

In this section, we will prove Theorem 12.11 Here, we will verify weak convergence via 
the convergence of finite-dimensional distributions and tightness. 
We first verify the convergence of finite-dimensional distributions. 

Lemma 3.1 The family of stochastic processes {X n (t), t € [0, 1]} converges in the sense 
of finite- dimensional distributions to the process X H defined by \2.J$. 



Proof: In order to prove Lemma l3.1j, we first establish some facts. 

Let us consider an arbitrary sequence of partitions of the interval [0, 1] of the form 



7T : U = to < ti < t 2 < ■ ■ ■ < t m = 1, 

with |7r m | — > 0, uniformly as m — > oo. Define 

in 

X m (t) = J2 K H(t,tT- 1 )W(A i ), (3.1) 

k=l 

where A< = [t^L v i™) and 

W(Ai) = W(tf) - W^). 

Fact 1. 

X m {t)^X H (t), (3.2) 



L 



0, where =^ denotes convergence in L 2 {Q). 
Since Kn(t, s) is continuous and monotone in s for every t, we can easily get that Fact 
1 holds. 



Define: 



where 
Fact 2. 



X m , n (t) = K n(t, C 1 )W n (A i ), (3.3) 

i=l 

W n (A i ) = W n (t i )-W n (t i - 1 ). 

X m , n {t) % X m (t), (3.4) 



as n — > oo, where denotes weak convergence. 

From the invariance principle and the continuous mapping theorem (see e.g., Billingsley 
[2] ), one can easily get that Fact 2 holds. 

Set 

X n {t) = } n K H {t,u)du—=; 

1 = 1 n 

Y nt \ r k m £ 

x m ,„(t) =y> I vxH(t,c 1 )i Ai(u) du-^=. 

k=i J ~^r i=i v 
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Furthermore, set 



K m (t,s) = J2KH(t,t™ l )lAM- 



(3.5) 



i=i 



Then if m (t, •) converges to K(t, •) in L 2 [0, 1]. 
Fact 3. As m — > oo, 

X m ,n{t) => X n (t). 



(3.6) 



Indeed, by E[£i£j] = if i ^ j, E[£f] = 1 and the Cauchy-Schwartz inequality, we obtain 



E 



X m ,n{~k) — X n (t) 



[nt\ 

Ef^fn / '" ^ (K m (t,u) - K H {t,u))du 

k=l 

[ntj 



& I 2 



r? 



k=l 
[ntj k 



fc-i 



k=l 
-1 



^E L (x m (*,«)--K - H(t,u))' 



du 



asm-) oo. 
Fact 4. For fixed n and t 



< / (K m (t,u) - K H {t,u)) du^O, 
Jo 



i? 



X m ,n(t) — X mtn (t) => 0, 



(3.7) 



as m — >■ oo. 

Let A is the Lebesque measure. Noting that J^i l^ fc (u)du = A([^-, ^) n A&), we can 

n 

rewrite X m n as follows, 



\nt\ 



t=l 



fc=l 



By (|3.3p . we can get 



A; — 1 fc N . \ ft 

,- nAiPi. 

n n / vn 



L"'/J 



i m ,n(t) = E^(Mr-i) E -%• 

1=1 fc=[nii-ij+l v 



(3.8) 



We first assume that 

Assumption: There exist {ii,j2"" >^m} ^ {lj **• > m } an d fo, • • ■ 
{1, 2 • • • , [nt\} such that ^) does not belong to Aj. and 



(3.9) 



n n 
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Noting that if C A,, A([^±, |) n A,) = I we get 



E 



Xm : n(t) X mn (t) 



E 



E 



o=l i=l 



T M 



1 = 1 0=1 * 



T 

E 

i=i 



n 



n q=l 



(3.10) 



From the assumption, it is easy to find a partition of [0, 1] with (£j ) g =i I ...,Af ) i-e., 

= *o < tji < tj 2 < ■ ■ ■ < tj u < tj M+1 = tj M+2 ■ ■ ■ = tj m = 1, 

and \tj q — tj J < |tj — for some i. 
Then, we get that as m — >• 00, 

Since Kn(t, •) is continuous and monotonic, we have that as m — > 00, 



E 



Xm t n(t) — X m ^ n {t) < / _ i 



i=l 



M 



(K H (t,j q ) - K H (t,j q ^)j (u)d 



8=1 



<C/ V^(t,j )lA dn^O. (3.11) 



By (|3.1ip . we get that the convergence is uniformly in n. 

On the other hand, if Assumption does not hold, then one can easily get that 

X m ,nify — X mn (^t). 

So Fact 4 holds. 

Fact 5. The family of processes {X n } converges to the process X H in the sense of finite- 
dimensional distributions. 

In order to prove Fact 5, it suffices to show that for any t±, ■ ■ ■ ,t p G [0, 1] and x € R, 
we have 



E 



exp(ix^2X n (ti)) ->• E exp (ixY X H (U)) 



(3.12) 



i=l 



i=l 



as n — > 00. 
Note that 



E 



exp(ix^2X n (ti)) — E exp (ra^J ff (t;)) 



i=l 



i=l 



<h + I 2 + h, (3.13) 
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where 



h 
h 



E 



E 



E 



exp [ix ^ X m , n (ti)) - exp (ix ^ X m (U) 

i=l i=l 
V P 

exp (ix ^2 X m>n (ti) \ - exp [ix ^ X n (t 



i=l 
P 



i=l 



ix^2x m (ti)J - exp (ix^2x H (tj 

i=i i=i 



We first deal with li- By the proof of Facts 3 and 4, we have that for fixed e > 0, the 
terms I2 can be bounded (uniformly in n) by e when m is large enough. Similarly, we can 
get that the same result for I3. For I±, from Fact 2, we can get that it goes to zero as 
n — > 00. From above arguments, we get that f|3. 12|) holds. 

Fact 6. 



E 



X n (t) - X n (t) 



0. 



as n — > 00. 
Indeed, 

E 



X n {t)-X n {t) 



/ fc-i n v n 



fc=i 

[ntj r 



k=l 



I nt I 

fc-i n 



< y (^(t, u) - K H (^,u)fdu ->• 0, 



(3.14) 



since Kn(t, •) is a continuous function. 
Next we will use Facts 1 to 6 to prove Lemma 13. 11 

In order to prove Lemma 13. 1\ it suffices to prove that for any t\ , • ■ ■ , t p € [0, 1] and 
rj G 1, we have 



E 



ir]y^X n (tj 

i=i 



E 



exp(ir)J2x H (ti) 



i=l 



(3.15) 



as n — > 00. 
On the other hand, we have 



E 



p p 
iT]y~)X n (ti)) - exp (ir]^X H (ti 



i=l 



i=l 



< D1 + D2 + L> 3 + £>4 + A>, (3-16) 
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where 



Di =E 



Do=E 



Do =E 



Da = E 



p 

i=l 
P 



exp ir] 



i=l 



exp (in ^ x m,n{ti) \ - exp (in ^ X n {ti) 

i=l i=l 
P P 

exp (irj ^2 X m ,n(tij) ~ exp (in ^ X m>n {t 



i=l 

p 



p 



exp (iv^2 X m>n (ti)) - exp (iy^2 X m {U) 
i=l i=l 



and 



D, =E 



exp X m (ti)j -exp(iri^2x H (ti) 



i=l 



i=l 



Then, using the same method as the proof of Fact 4, we can prove (|3. 15j) holds. 



□ 



Next, we prove the tightness of {X n (t)} n£ jq. 
Lemma 3.2 The family {X n (t);t £ [0, 1]} given by \2.1^ is tight. 
Proof: In order to prove this lemma, it suffices to prove that for any s < t < r G [0, 1], 

\X n (t) - X n {s)\\X n {r) - X n (t)\] < C\r - s\ 2H . 



E 



(3.17) 



since H > ~. 

Noting that the kernel Kn(t,u) vanishes when u is larger than t, we have 



E 



X n (t)-X n {s) 



E 



\nt\ 

£ 

i=l 



n 



[nt\ 



[ns\ 



K H ( ,u) - K H ( ,u))du 



n 



n 



[nt\ 
i=l 



K H (^,u)-K H (^,u))du)\ (3.18) 
n n I / 



since E[&&] = if i ^ j, and E[f ?] = 1. 
By the Holder inequality, we get that (|3. 18)) can be bounded by 



LntJ 



i=l 



*-~f.li=l V n n Iq \ n n J 

[nt\ [ns\ 2H 



< C 



n n 



since (H^j) holds. 

Therefore, the tightness can be obtained by Billingsley Section 5]. 



(3.19) 



□ 
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Proof of Theorem \2.1\ Theorem 12.11 is a direct consequence of Lemmas 13.11 and 13.21 
because tightness plus convergence of finite dimensional distributions implies the weak 
convergence in the skorohod topology (see Billingsly [2j Section 5]). □ 
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